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Abstract In the fireball model, it is more physically realistic that gamma-ray burst (GRB) 
ejecta have a range of bulk Lorentz factors (assuming M cx F^"). The low Lorentz factor 
part of the ejecta will catch up with the high Lorentz factor part when the latter is decelerated 
by the surrounding medium to a comparable Lorentz factor Such a process will develop a 
long-lasting weak reverse shock until the whole ejecta are shocked. Meanwhile, the forward 
shocked materials are gradually supplied with energy from the ejecta that are catching-up, 
and thus the temporal decay of the forward shock emission will be slower than that without 
an energy supply. However, the reverse shock may be strong. Here, we extend the standard 
reverse-forward shock model to the case of radially nonuniform ejecta. We show that this 
process can be classified into two cases: the thick shell case and the thin shell case. In the thin 
shell case, the reverse shock is weak and the temporal scaling law of the afterglow is the same 
as that in Sari & Meszaros (2000). However, in the thick shell case, the reverse shock is strong 
and thus its emission dominates the afterglow in the high energy band. Our results also show 
slower decaying behavior of the afterglow due to the energy supply by low Lorentz factor 
materials, which may help the understanding of the plateau observed in the early optical and 
X-ray afterglows. 
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1 INTRODUCTION 

The central engine and surrounding environment provide the most important insights to the mystery of 
gamma-ray bursts (GRBs), the most violent explosions in the universe. Thanks to BeppoSAX, owing to its 
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ability to accurately locate objects, the first afterglow of a GRB was discovered in 1997 (Costa et al. 1997). 
Afterward, broad band data of afterglows were achieved, and were fitted using the standard fireball-shock 
model (Rees & Meszaros 1992, 1994; Meszaros & Rees 1993, 1997; for reviews, see Zhang 2007). The 
parameters of GRBs, such as the total burst energy, the type (ISM or wind) and number density of the en- 
vironment, and the electron and magnetic field equipartition factors were then constrained (e.g., Wu et al. 
2003; Fan et al. 2002; Zhang et al. 2003). Early afterglows can even be used to constrain the initial Lorentz 
factors of GRB fireballs (MoUnari et al. 2007; Xue et al. 2009). It is believed that early afterglows are 
produced by reverse-forward shocks when relativistic ejecta interact with the circum-burst medium, which 
was first studied by Rees & Meszaros (1992) and Sari & Piran (1995). Then GRB 990123, a remarkable 
event with a bright early optical flash was discovered (Akerlof et al. 1999), which was interpreted by the 
reverse-forward external shock model weU(Sari & Piran 1999). Much attention was consequently focused 
on early afterglow radiation, considering the effects of circum-burst environments and non-relativistic re- 
verse shocks (Kobayashi 2000; Wu et al. 2003; Zou et al. 2005). Subsequently, more optical flashes were 
observed in, e.g., GRBs 021211, 050525a, 060111B, 060117B and 080319B (Wei 2003; Shao & Dai 2005; 
Klotz et al. 2006; Jelinek et al. 2006; Racusin et al. 2008). These works were based on the assumption that 
the Lorentz factor does not change in the shell. The whole Ught curve from the reverse-forward external 
shock has two types: re-brightening (Type I) and flattening (Type II) (Zhang et al. 2003). In some GRBs, 
such as GRBs 050319, 060206, 060210 and 060313, the early optical light curves have a plateau which is 
difficult explain within the uniform ejecta model (UEM). Meanwhile, the shaUow decay of the canonical 
X-ray afterglow behavior discovered in the Swift era (Nousek 2006; Zhang et al. 2006; O'Brien et al. 2006) 
remains a matter of debate. 

Because of the above problems, we reconsider the bary on-dominated energy injection model in which 
the ejecta have a wide F-distribution: the part of the ejecta with low Lorentz factor lagging behind the one 
with high Lorentz factor. The low-F part catches up with the high-F part when the latter is decelerated 
to a comparable Lorentz factor, so the reverse shock is usuaUy mildly relativistic and mainly contributes 
to the far-IR or miUimeter band (Rees & Meszaros 1998; Sari Meszaros 2000). Once the reverse shock 
starts, it wiU travel through the whole ejecta from the front highest-F part to the rear lowest-F part. Based 
on the treatment widely adopted in the UEM, we reconsider this issue by assuming a given distribution of 
Lorentz factors in the ejecta. We calculate the dynamic evolution of the reverse-forward shocks produced 
by this radially structured ejecta propagating into the circum-burst medium, and present the analytical and 
numerical results. 

We organize our paper as foUows. In §2, we describe the dynamics of the reverse-forward shocks in- 
cluding the thick shell case and the thin sheU case, respectively. We discuss the reverse shock emission in 
§3. The numerical results are shown in §4. FinaUy, we present a brief discussion in §5. 

2 DYNAMICS OF THE EJECTA WITH A F-DISTRIBUTION 

As Rees & Meszaros (1998) postulated, the central engine of GRBs may eject relativistic shell-hke ejecta 
with a range of Lorentz factors 
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Such ejecta have an energy TM oc r^*+^. Applying the model to observations shows that the value of the 
index s is typically ~ 2.5 (Zhang et al. 2006; Nousek et al. 2006), which is larger than the suggested value 
^ 1.5 (e.g., Rees & Meszaros 1992), so the low-F mass carries more kinetic energy than predicted. When 
the ejecta interact with the circum-burst medium, a pair of shocks emerges: a forward shock propagating into 
the circum-burst medium and a reverse shock propagating into the shell. There are four regions separated 
by the two shocks: (1) the unshocked circum-burst medium, (2) the shocked medium, (3) the shocked 
shell material and (4) the unshocked shell material. Using the shock jump condition and the equality of 
pressure and velocity along the contact discontinuity, the Lorentz factor F and the number density in both 
shocked media can be determined by the density of the circum-burst medium ni and the unshocked shell 77,4 
(Blandford & McKee 1976, hereafter BM). Here, the number density of the unshocked shell is nonuniform, 
depending on the Lorentz factor distribution in the shell 

^ \{dM{>r)/dT)idr/dx)\ ^ 



We assume an initial Lorentz factor distribution in the shell 
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where Aq and T^i^ are the initial width and minimum Lorentz factor of the shell, x{t — 0, F) represents 
the initial position in the shell with the origin located at the outer edge of the shell (See Fig.[T]i- Due to the 
distribution of Lorentz factors, the shell will spread with time, then the value of x of a fixed element will 
become larger and larger, 

x{t, F) = x{t = 0, F) + (/3„,ax - f3)ct, (4) 



where (3 — \J\ — l/F^. Thus the Lorentz factor distribution in the shell at any time is determined. 
The properties of the shocks are largely determined by the parameter defined as 
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where Mi, is the total mass of the shell. Eq. (1) can be written as A/(> F) « Af^F^jj^F^* (here Fmin ^ 
Tmax is assumed). Combining Eq. (HJl and (|5]l, we need another equation to describe the evolution of the 
radius and Lorentz factor of the shocks: the relation between the distance dx which the reverse shock travels 
in the shell and the distance dr which the shell propagates in the circum-burst medium in the same time 
interval is (Kobayashi 2000) 

dr = aVf^dx, (6) 

where r is the radius of the shell and the parameter a is '-^ 1. 

There are two approximations under which the shock evolution can be described analytically: the thick 
shell case and the thin shell case, depending on the significance of the spreading effect in the last term in Eq. 
(|4]l. If the spreading term is larger than the initial x{t — 0, F), the shell is regarded as a thin shell, otherwise 
it is a thick shell. Below, we will consider these two cases separately. 

2.1 Thick Shell Case 



In the thick shell case, the spreading effect can be ignored, so the width of the shell always remains at its 
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We can get from Eqs. (O and ^ 



dx = — 6Ao 



-Mr 



(7) 



We can now calculate the comoving number density 714 and the density ratio /. According to Eq. (|5]l, we 
have 

dr ^ ahM,Ti^,^T-^ dT _ 
dx A'Kr'^mpTii dr 

In general, the number density of the circum-burst medium can be modeled as ni — Ar^*' . Specifically, 
A = ni = 1.0 riifi cm^^ for an ISM environment (fc = 0) and A = 3 x 10"^^ A^, cm^^ for a free wind 
environment (fc = 2) (Chevalier & Li 2000). Throughout this work, we adopt the convention = Q/IQ^ 
in cgs units. The solution of Eq. ^ reads 
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Keep in mind that dr — 2T^cdt, the evolution of the Lorentz factor and radius of the ejecta with time can 
now be described. Just when the reverse shock crosses the shell, the Lorentz factor is equal to Fmin and the 
shell reaches the crossing radius 
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In the thick shell case, the spreading effect is always negligible which requires ta < Fj^jj^ Aq. Thus a lower 
limit of the minimal Lorentz factor must be satisfied 
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Whether or not the reverse shock is relativistic depends on the parameter 
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It can be seen that the parameter //F^ is smaller than one for the typical wind case (k=2) and for a 
dense ISM case (rii = 100) at the crossing time, and before the crossing time it is proportional to 
p2(3b-6fe+2fe-s-7)/(4-fc) jjjg Lorcntz factor distribution is not too steep < & < 3 and the mass dis- 
tribution index has a typical value of 1.5 < s < 2.5, f < holds for the entire reverse-forward shocks 
interaction period which means that the reverse shock is relativistic all along. While in the other parameter 
space, it is possible that / 3> F^ at the initial stage when the reverse shock is non-relativistic, and then 
/ < F'^ which corresponds to the reverse shock evolving from being non-relativistic to being relativistic. If 
the reverse shock is relativistic, the relative Lorentz factor 734 between the shocked shell and the unshocked 



Early afterglows from radially structured outflows and the application to X-ray shallow decays 



5 



2.2 Thin SheU Case 

In the thin shell case, the spreading effect is dominant. The position of an element in the shell, x, can be 
approximated by 

x{T) « (/3„,ax - (i)ct « (13) 

the parameter / can be written as 

is-i)E ( r 

Based on the same procedure appUed in the thick shell case, we can obtain the relation between the radius 
and the Lorentz factor 



Then the crossing radius is 

'3-fcN2 o?{s-\)E 



rA = 



47rAmpc2r2.„ 




'^'EU\l},/'rJi%cm, (fc = 0) 
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E53A-^T-l^^ cm, (fc = 2) 

In contrast to the thick shell case, an upper limit of the minimal Lorentz factor must be satisfied to keep the 
thin shell assumption valid all along 

■ Ao:?fn-;/«, {k = 0) 

1/4 ■ (17) 

-^53 ^0,12 ^* > {k = 2) 

For the thin shell case, it is interesting that //F^ = a~^(s + 1/3 — fc)^, which means that the reverse shock 
is always mildly relativistic (734 — 1 w F^// ~ 1). 

3 REVERSE SHOCK EMISSION 

Now that the dynamic related parameters, i.e., F, 734 and 714 are determined, the radiation related properties 
of the shocked materials, such as the strength of the magnetic field S ■ , the minimum Lorentz factor 7^^^, the 
cooling Lorentz factor 7^ ^ and the number N^^i of shocked electrons can be determined. For the shocked 
region, the fraction e b and Ce of the internal energy are assumed to be carried by magnetic fields and shock- 
accelerated electrons, respectively. The co-moving magnetic field is equal to -^/STrese^, where the internal 
energy density = (7rei — l)mj,c^n, (for the forward shock 7rei = 72 ~ 1/2F^/^/^/*; for the reverse 
shock 7rei = 734). The minimum Lorentz factor 7^ • = eeCp(mp/me)(7rei — 1) withCp = (p— 2)/(p— 1), 
the cooling Lorentz factor 7^ ^ = dirniec/ (arBf-yit). The increase of the number of shocked electrons is 
dNe,3 = 4:nr^rn4dx for the reverse shocked region and dNe^2 = 47rr^nidr for the forward shocked 
region. In the standard synchrotron radiation model, the two characteristic frequencies and the peak flux 
density are 

QeB'i ,2 _ ge-B- ,2 TP _ NeJ_meC^aT 



o -7m,i7i, i^c,i = 7^ 7c,i7i, -Pl^.max.i = ■Bi7i, (18) 



where dL is the luminosity distance of a GRB, qe is the charge of electron and ar is the Thomson cross 
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listed in Table 1 for both the forward shock and reverse shock, and for both the thick shell case and thin 
shell case. 

The distinct discrepancy between the nonuniform ejecta model (NUEM) and the UEM is the reverse 
shock emission. Once the reverse shock has crossed the shell, the forward shock and shocked region begin 
to approach the Blandford-McKee (BM) solution (Kobayashi et al. 1999), and the following light curve is 
the same as that in the UEM. Below we only discuss the synchrotron emission from the reverse-shocked 
region before the crossing time. Since we have already obtained the temporal indices of i/„j, Vc and F^ax, 
in the following we only need to know the values of the characteristic frequencies and the peak flux density 
at the crossing time so we can extrapolate the early light curve back in time from the reverse shock. 

3.1 Thick Shell Case 

The reverse shock in the thick shell case we consider here is assumed to be relativistic (734 — 1 « 734). The 
crossing time Ta is ~ A/c, when the Lorentz factor of the shell is Fmin and the number of the shocked 
electrons is the total number of electrons in the ejecta, i.e., iVe_3 = E{s — 1)/ {smpC^T^i^). According to 
Eq. (fTSl i. we have 

v^^A.Q X 10^3(1 + z)-i(^)'e^ _o.56l/i2n;grLn,2.5 Hz, (19) 



1.0 X 10^^(1 + 
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for the ISM case, and 
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for the wind case, where Ta,2 = T/lOOs. The above expressions for the synchrotron radiation at the 
crossing time are quite similar to those in the UEM (e.g., Kobayashi 2000; Wu et al. 2003). We find that for 
a set of combinations of reasonable parameter values (s ~ 2.5, 2, 1.5 and b = 2,1) the tcm (the time when 
i^m = Vc) is always small for the ISM case, indicating that the reverse-shocked electrons are always in the 
slow cooling region. However, for the wind case, the electrons are usually fast cooling during the entire 
reverse shock phase, because Vm is typically larger than Vc at the crossing time and Vm decreases with time 
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3.2 Thin SheU Case 

In the thin shell case, the reverse shock is always inildly-relativistic (734 — 1^1) and the crossing time 
depends on the crossing radius as t/^ oc r^^'^^~'^''^^^^~^^\ For convenience, we choose s = 2 to give the 
typical values of the two characteristic frequencies and peak flux density of synchrotron radiation, 



2.1 X 1012(1 + ^)-i(^)Ve^__o.5e^/i2rii„_i.8 Hz, (25) 

^ 7.1 X 101*^(1 + z)-ia-^/^e^;'Z',i;3-//^n-^/«r^t,,.3 Hz, (26) 

F.,ma^ ~ 2.3(1 + ^)£>2"8'eB,-2-S53n}fo r^in.l.S Jy, (27) 

for the ISM case, and 

i^m - 4.2 X 10^2(1 + 2)-i(^)'a2e2__o.5e^/i2A2f!o.5i55-3irL„,i.8 Hz, (28) 

i/, ~ 1.9 X 10^3(1 + ^)-'a'e-'/|E53^:,!/o.5r-?n,i.8 Hz, (29) 

F,,max - 385.0(1 + Z)i?2t«"'4(-2^!f-0.5rmi„,1.8 Jj. (30) 



for the wind case. It is shown in Table 1 that the temporal indices of Um, Vc and Fmax are the same as in the 
varying injection model (Sari & Meszaros 2000), which indicates that the F-distribution in the shell does 
not affect the shape of the Ught curve for the thin shell case (see Table 2) because the spreading effect erases 
the initial F-distribution. 

4 NUMERICAL RESULTS 

Since the thin shell case in our treatment has the same results as that in Sari & Meszaros (2000), hereafter 
we just present the numerical results of the thick shell case. For the sake of simplicity, we only show the 
result by one set of parameters with s = 2 and 6 = 2, which is nevertheless sufficient to compare the result 
of the thick shell case in our paper to that of the thin shell case in the literature. We assume that the redshift 
of a GRB is z = 2 since it is the average value of the observed GRBs in the Swift era (Le & Dermer 
2006). According to the standard shock acceleration mechanism, the energy index p of electrons is about 
2.2 ^ 2.3, so we choose p = 2.3. 

We follow the method of Zou et al. (2005) to perform our numerical calculations with Fmax = 1000, 
Fmin = 300, Eo = 1.0 X 10^^ erg, Aq = 1.0 x lO^^ cm, Ce = 0.3, cb = 0.01, A = 0.1, m = 10 cm'^ 
and adopt the standard cosmology model with Qm = 0.27, CIa = 0.73 and Ho = 71 km Mpc"^. 

Figs. 2 and 3 show the light curves of synchrotron emission in the optical R band and X-ray band. To 
compare with the results of the UEM, we also plot the light curves of the UEM in the same energy band. 
In our calculations, the nonuniform and uniform ejecta have the same kinetic energy and total mass, so 
the uniform ejecta have a Lorentz factor of sFmin/ (s — 1) = 600. The upper panel denotes the R-band 
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(1) In these two energy bands, there is a slow decay phase before the crossing time for the NUEM which 
is attributed to the energy injection by the low-F part. (2) Before the crossing time, the reverse shock 
emission dominates the radiation in the optical band, while in the X-ray band, the shocks have comparable 
contribution. After that the forward shock emission gradually becomes important. This result is similar to 
the UEM. It is possible that the two shocked regions may have different microphysical parameters and 
es, which do not change the first conclusion but might change the second one. 

Although we do not show the Ught curves of the thin shell case, we would like to emphasize the dif- 
ferences between the thick and thin shell cases. Since the reverse shock in the thin shell case is mildly- 
relativistic all along and the flux density has the same behavior as presented by the refreshed shock scenario 
(Sari & Meszaros 2000), a simple test of the thin sheU case is that it predicts a maximal flux in the far- 
infrared or milUmeter range a few hours to a few days after the GRB trigger. The forward shock emission 
contributes mainly in the high energy band and decays slower than that of the normal standard forward 
shock model. However, in the thick shell case, the reverse shock may be relativistic and dominate the flux 
in both the optical band and X-ray band at early times. 

5 DISCUSSION 

We have described the dynamics of radially structured ejecta interacting with the circum-burst medium 
through extending the method used in the UEM. It can be classified into two types: the thick shell case 
and the thin shell case, of which the latter is the same as in the UEM. Two parameters are introduced in 
the NUEM. One is the mass distribution index s and the other is the Lorentz factor distribution index b. In 
our treatment, the thin shell case reproduces the same results obtained by Sari & Meszaros (2000). On the 
other hand, in the thick shell case, the reverse shock could be relativistic or initially be non-relativistic and 
then become relativistic, contributing comparable radiation in the high energy band as the forward shock. 
Anyway, the energy injection induced by the nonuniform ejecta sweeping up the surrounding medium 
causes the light curves to decay more slowly. 

Observationally, GRB990123 was seen to have a bright optical flash with initial flux decay as F oc t~^, 
which is attributed to the reverse shock emission, and subsequently the optical afterglow decays as oc 
t~^'^, which is mainly due to the forward shock emission. Both the thick and thin shell scenarios in the 
UEM can fit this optical flash well. The fast rise of t"^-^ can be explained if the circum-burst environment is 
homogenous ISM (Kobayashi 2000; Fan et al. 2002). However, some optical flashes (e.g., GRBs 021211, 
050525a, 0601 1 IB and 0601 17B) have not been observed with this early rising part. This may be intrinsic, 
or due to late responses and slow slewing of optical telescopes. An early optical plateau was observed in 
a few GRB afterglows, e.g., GRBs 050319, 060206, 060210 and 060313, which may be attributed to flie 
relativistic reverse shock emission of the radially structured ejecta sweeping up the circum-burst medium. 

In Swift GRB afterglows, pecuUar chromatic breaks have been universally observed, but the origin of 
these chromatic breaks is stiU an open question. Some models were proposed to explain the chromatic 
breaks, such as the model with evolving microphysical parameters, or the model in which the optical and 
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involving reverse shock emission (Genet et al. 2007; Uhm & Beloborodov 2007). Our current work shows 
that the chromatic breaks cannot be due to the reverse shock and we favor the former explanations. 

Although the detailed prescription of a radially structured shell propagating into the circum-burst 
medium is presented in this paper, it should be noted that the shock-heated material separated by the contact 
discontinuity is assumed to be uniform. A more accurate solution of the reverse-forward shock interaction 
and emission needs a relativistic hydrodynamic simulation. 
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2(7-2/c+s) 


1h{A-k) 


2(7-2/c+s) 



Table 1 Temporal indices of the peak frequency v^, the cooling frequency v^. and the peak flux 
density -F,^,max for both the forward shock and reverse shock. 



thick shell case 




V < min[t'ni,t'c] 


min[i/m,t'c] < f < maxfi/nnt'c] 


max[fni,i'c] < f 




Slow Cooling 
Fast Cooling 


5-96-fe+26fe+15s-4fcs 


12-126-3/c+36fe-16p+5fep-6fep+20s-5fes-fcps 


16-166-6fc+66fc-16p+5fep-6fcp+16s-2fcs 


— kps 


36(4-fe) 
-5-76+3fc+17s-6fes 


46{4-fc) 
-16b-k+5bk+16s-3ks 


46{4-fc) 

16-16i)-6fc+6i)fc-16p+5fep-i)fep+16s-2fes 


— kps 


3b(4-k) 


4b(4-fc) 


46(4- A:) 




thin shell case 


Slow Cooling 
Fast Cooling 


12-3k-9s+4ks 


12-3fc+12p-3fcp-123+6fc3+*:ps 


16-6fc+12p-3fcp-8s+2fcs+fcps 




6fc-3{7+s) 
16-3fc-lls+6fcs 


8fc-4{7+s) 
28-9fc-8s+3fes 


8fe-4(7+3) 
16-6fe+12p-3fcp-8s+2fcs+fcps 




6fc-3{7+s) 


8fc-4{7+s) 


8fc-4(7+s) 





Table 2 Temporal indices of the flux density F^, oc t^u^ of synchrotron radiation from a reverse 
shock. Both the thick sheU case and the thin sheU case are considered. 




O 

Fig. 1 Schematic of the radially structured shell. A is the width of the shell, Fmin and Fmax are 
the Lorentz factors on the edge of the shell respectively. The part with the lower Lorentz factor 
in the sheU has more kinetic energy (marked with darker gray) than that with the higher Lorentz 
factor. 
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Fig. 2 Synchrotron radiation flux density in the R band and X-ray band as a function of time 
for the ISM case. Solid lines represent the whole emission from both the forward shock and the 
reverse shock. Dashed lines represent the contribution from the reverse shock. Thick lines are for 
the NUEM while thin lines are for the UEM. 



Xue-Wen Liu, X. F. Wu, Y. C. Zou & T. Lu 




Fig. 3 Synchrotron radiation flux density in the R band and X-ray band as a function of time, the 
same as in Fig. 2. for the wind case. Parameters are given in the text. 



